Directional (version 2.2)

Mle of the Matrix Fisher distribution on SO(3): Mle of the Matrix Fisher distribution on SO(3)

Description

It returns the maximum likelihood estimate of the Matrix Fisher parameter F(3x3).

Usage

matrixfisher.mle(X)

Arguments

X

An array containing rotation matrices in SO(3).

Value

The components of $ svd( \bar(X) ) $.

References

Prentice M. J. (1986). Orientation statistics without parametric assumptions.Journal of the Royal Statistical Society. Series B: Methodological 48(2).

See Also

rmatrixfisher

Examples

F <- 10^(-1) * matrix( c(85, 11, 41, 78, 39, 60, 43, 64, 48), ncol = 3 )   ### An arbitrary F matrix
X <- rmatrixfisher(3000, F)
matrixfisher.mle(X)
svd(F)
## Bigger sample is needed to obtain better results
## X <- rmatrixfisher(30000, F)
## matrixfisher.mle(X)
## svd(F)